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We used quantum process tomography to investigate and identify the function of a nonideal 
two-qubit quantum-state filters subject to various degree of decoherence. We present a simple 
decoherence model that explains the experimental results and point out that a beamsplitter followed 
by a post-selection process is not, as commonly believed, a singlet-state filter. In the ideal case it is 
a triplet-state filter. 
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PACS numbers: 03.65.Wj, 42.50.Dv 

Characterization of quantum devices is important, and 
the characterization of single- and two-qubit devices is 
particularly important because a single-qubit rotation 
gate and a two-qubit controlled-NOT (C-NOT) gate are 
two building blocks of the quantum computer|lj. The 
quantum systems in actual devices may undergo nonuni- 
tary state evolution because they interact with their en- 
vironments, so identifying and understanding the func- 
tion of such nonideal devices are of practical impor- 
tance. Nonideal devices can be characterized by quan- 
tum process tomography (QPT)0,|3,0], and this method 
has been used to characterize single-qubit transmission 
channels^, and rotation gates for photonsQ, gates 
for ensembles of two-qubit systems in NMRj8j, and a 
two-qubit quantum-state filter and a C-NOT gate for 
photons |Toj. 

This letter reports the QPT of a quantum-state filter — 
the beamsplitter (BS) in a Hong-Ou-Mandel (HOM) 
interferometer ^3, ^2 — as a demonstration how the func- 
tion of two-qubit devices can be identified experimen- 
tally. The result of an earlier investigation of the nearly 
ideal filter @ was inconsistent with a common belief that 
an ideal BS followed by a post-selection process — only 
the outcomes where the two input photons come out in 
two output ports are taken into account — behaves as a 
singlet-state filter 0, but the reason for this in- 
consistency was not fully understood. To resolve this 
problem, we have investigated the function of the non- 
ideal filters subject to various degree of decoherence 
and have identified two leading contributions to the 
operator-sum rep resentation of the quantum operation 
for them0,H . Here we present a model that explains 
our experimental results and resolve the above inconsis- 
tency. 

Let us start by introducing the process matrix |To). 
which offers us the most general description of the quan- 
tum device allowed in quantum mechanics, and showing 
how to obtain it experimentally. Consider a two-qubit 
device whose function is to map a two-qubit input den- 
sity operator pi n to a two-qubit output density operator 
Pout = £{f>in)i where £ is a completely positive linear 
mapQ. The map £ is considered a superoperator act- 
ing on positive operators Let us introduce the set of 



standard operator basis {Xk} k=0 that spans two-qubit 
operator space^], where Xuj\ = = 0,1) and 

(ij) := 2i + j. Then £ can be expanded as 



m = 
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where [kl] := 4k + Z |T5| . A 16x16 positive Hermitian 
matrix x S fully characterizes the function of quantum 
device and is called the process matrix in the standard 
operator basis jicj. 

We can transform \ m t° the- process matrix in any 
operator basis set {^al^Lo t na f spans the two-qubit op- 
erator space. Note that A a may be defined either locally 
or nonlocally. For example, the Kronecker products of 
the Pauli operator basis B[m] = Ej~®Ei (k, I = 0, • • • ,3), 

where {E k }={a Q / V2, o\j\[2, &2/V2, a 3 /V^} and a de- 
notes an identity operator, is a local operator basis, 
whereas the outer product C[ki] = |$fc) where |$fc) 
are Bell states, is a nonlocal one. The operator basis A a 
is, in general, unitarily related to the standard operator 
basis X k <g> X t . That is, A a = Y,l,l=o {U) [kl]a X k ®X h 
where U is a 16x16 unitary matrix. Then £ can be ex- 
panded using A a as follows: 



m = E U A ) 

q,/3=0 



a/3 AaPAp, 



(2) 



where \ A — U^x S U . Therefore, \ A is given by the uni- 
tary transform of x S ■ 

To obtain \ s , it is sufficient to determine the as- 
sociated four-qubit state 01- Let us consider a four- 
qubit system composed of the two two-qubit systems 
1-2 and 3-4 and introduce the double-ket notation 
|7)) := (7 <g> 7) |$) to denote the bipartite state, where 
|$) = nH,i) is the unnormalized maximally en- 
tangled state H, [lg • Let us introduce the unnormal- 
ized four-qubit density operator T>£ associated with £ 
by V £ = £^ ®Z( 24 )(|7)) 12 «7| ® |7»34«7|), where 1 
is an identity map and the superscript and subscript re- 
spectively denote the space on which each map acts and 
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the space to which the operator belongs |l7j. Using Eq. 
^and \Xk)) := (Xk ® -01$) j we can rewrite £>g as 



E (X 5 ) [fei][mn] l^»l^»((^l((^|. (3) 



k,Lm,n— 



If we note that |-X7#\)) = Eq. ^indicates that \ S 

agrees exactly with the density matrix of the associated 
state T>£ given in the standard state basis. Evaluation of 
the process matrix is thus reduced to evaluation of the 
associated state T>£. 

To obtain X>g, it is sufficient to evaluate all the 
maps £{Xuj\ ® X/)-l)) for all the elements of the stan- 
dard operator basis X^ <g> X^i), since we can ob- 
tain the permutation equivalent of T>£ — that is, T>g = 

P23P12PUV £ P M P 12 P 2Z = X( 12 ) ® f( 34 )(|/)) 13 ((/| ® 

\I))2a{{!\) — by using the following identity: 



<Xi®£(X k ®Xi), 



(4) 



fej=0 



where Py is a permutation operator between the sys- 
tems i and j. To evaluate these maps, we need 
only to use state tomography^^ to evaluate the 16 

two-qubit output states {^{p'p ® pi- )}|,j=o associated 
with the set of 16 separable input states, where pi £ 
{|0)(0|, |1)(1|, |+)(+|, |0)(0|}, |+) = l/v / 2(|0) + |l)), and 
|0> = 1/V2(|0> + i\l)) HH0. We therefore can obtain 
the process matrix directly from the experimental data 
without using a matrix inversion procedure Q, Q. 

We used the above procedure to obtain the pro- 
cess matrices for nonideal two-qubit quantum-state fil- 
ters. The experimental setup is shown in Fig. ^ 
A 0.3-mm-thick type-I crystal (BBO) was pumped by 
frequency-tripled pulses generated by 100-fs pulses from 
a mode-locked Tksapphire laser (repetition rate «82 
MHz). Pairs of horizontally polarized 532-nm photons 
were generated and those traveling in two directions 9° 
from the pump beam were selected by 5-mm irises (1 m 
from the crystal) and identical 3-nm-width filters (IF) 
placed before the detectors. They were injected into a 
HOM interferometer [ti ll ll Il2| having a broadband non- 
polarizing cube BS (Suruga Seiki model S322-20-550). 
The horizontally (H) and vertically ( V) polarized states 
(respectively |0) and |1)) and the other two states re- 
quired for the QPT were prepared by half- and quarter- 
wave plates (HWP and QWP) placed immediately before 
the BS. Down-converted photons were incident on the BS 
at 45°. After the BS, polarization analyzers consisting of 
a QWP and a HWP followed by a Glan-Thompson prism 
(POL) and a single-photon counting module (PMT) per- 
formed projective measurements onto the set of the 16 
product states {pp ® p*- }? J=0 . Only signals coincident 
at the two PMTs were recorded, and the associated 16 
output states \£{pi ®Pj )}f iJ=0 were evaluated by state 
tomography 0| . Note that we used unnormalized out- 



put density matrices because the state filtering is trace- 
decreasing process whose output intensity depends on the 
input state. 
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FIG. 1: Experimental setup to investigate the beamsplitter. 

We used optical trombones to equalize the two path 
lengths to within the coherence length l c of the down- 
converted photons (~ 200 p,m). Two two-photon am- 
plitudes contributing to coincidence detection — i.e., that 
both photons are reflected (r-r) or that both are trans- 
mitted (t-t) — were initially (t = 0) aligned to overlap 
both spatially and temporally, giving a HOM dip visibil- 
ity of 85% as shown in Fig. [3 In this case, the BS is 
commonly thought to work as a singlet-state filter [T3 . fl3| - 
To investigate the nonideal filter, we introduced a delay 
r in one arm by moving one of the trombone prisms. 
The larger the r is, the smaller the overlap of the two 
amplitudes is and the more the temporal and polariza- 
tion information of the photon pair is entangled. We will 
show later that this results in decoherence if we observe 
only the polarization information. 
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FIG. 2: HOM dip observed in the coincidence record. Arrows 
indicate delays r of (A) 0, (B) 100, and (C) 350 fs introduced 
in one of the two arms. 

We obtained the process matrices \ F associated with 
three-delay conditions (r = 0, 100, and 350 fs). To make 
the results clear, we represent them here in the nonlocal 
operator basis = (Ei ® Ej)Us, where Ug is a swap 
operator. They are shown in Fig. [31 where space is saved 
by showing only the real parts. What is most interest- 
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ing is that the contribution of (x F )isi5 increases with 
increasing delay r. 

This result is consistent not with the common belief 
but with what was pointed out by Mitchell et a/.|9(. In- 
deed, we found that our results are well reproduced if 
we assume that £(p) has an operator-sum representation 
with two Kraus operators V~ and 'P + (lHl4|: 

£{p) = (1 - p)fi-pf>-^ + P V+pV + \ (5) 

V ± =TUi±nU 3 (0 1 ,0 2 )Us, (6) 

where p (0 < p < 1/2) is a parameter that can 
be interpreted as the degree of decoherence and where 
U 3 (di,9 2 ) = e l8ia3 / 2 a 3 ® e- ie2 * 3 / 2 <7 3 , which explains 
the large contribution of (x F )i5 15 instead of (x F )oo 
as would be expected for the nonidcal singlet-state fil- 
ter. The Kraus operators V~ and V + correspond re- 
spectively to the operator for the ideal state-filter and 
the noise operator. Note that Kraus operators are not 
unique and are not necessarily orthogonal Ij. Figure 0] 
shows the process matrices calculated using the measured 
value 7Z/T = 0.76 and appropriate values of the param- 
eters 0i, 6*2, and p. As shown in Fig. Eqs. |S]and[f)] 
well reproduced the experimental results. It follows that 
the operator V~ describing ideal operation deviates from 
the singlet-state because of the factor IJ 3 (^\, Q%), which 
is consistent with Ref. . 

To understand these results, we need a decoherence 
model. Let us first call our attention to the function of 
the BS in the polarization space. Let a and b (c and d ) 
denote input (output) field operators of the BS, and the 
polarization be indicated by subscripts (e.g., by an and 
ay for the H- and V-polarized field operators). If we use 
a right-handed coordinate system defined by the positive 
directions of H-, V- polarized field (p and s components) 
and propagation vector (Fig. ^) to describe the input and 
output fields, their operator equations are given by: 

c H = Vfa„ + ie i( - r+ ^Vni> H , 

d H = ie-^+^VK&H + Vfb H , 

cy = VTav + ie lA vTZbv, 

d v = ie- lA Vnd v + VTbv, (7) 

where T and 1Z are respectively the transmission and re- 
flection coefficients of the BS. We have assumed that the 
p and s components of the field respectively acquire phase 
shifts r + 7r and A upon reflection. The factor ir accounts 
for abrupt change in p-s phase shift due to nonadiabatic 
changes in the k vector upon reflection, which has been 
frequently noted in the papers concerning the topological 
phase effect in optics [2(J, |2l|, |23 • Note that the helicity 
of a photon in a state |0) changes sign upon reflection 
from the BS but not upon transmission through the BS if 
A — T = 0. If we consider the output fields that can con- 
tribute to coincidence measurement, we conclude from 
Eqs. [7| that the BS projects the input state onto V~ 
with 6i = 9 2 = A — r, which may be finite for a non- 
polarizing BS. Note that the projector V~ is the coherent 



sum of the identity operator Ui associated with trans- 
mission and the operator U3(0i,0 2 )Us associated with 
reflection, where 1/3(0%, 62) accounts for the helicity re- 
versal and the extra p-s phase shift upon reflection. 

Although the above description may suffice when r = 
0, the temporal (T) space of the two-photon state should 
be taken into account when t ^ 0. Concerning the tem- 
poral state space, it would be reasonable to assume that 
the BS leaves the state untouched on transmission and 
swaps the state on reflection. We thus assume that the 
function of the BS is described by the following operator 
entangled with respect to the polarization and temporal 
state spaces: 

f PT = TUf g> Uj -7111^(9^02)11% <» L>J, (8) 

where the superscripts indicate the space on which each 
operator acts. We further assume that input state is un- 
cntangled, i.e., that it is a product state of polarization 
state p and temporal state w. Then, by partially trac- 
ing over the temporal degree of freedom, we obtain the 
following output polarization state £(p) consistent with 
Eqs. E] and HO 

£(p)= Tr T V PT (p®u) (P PT ) ] 

= (I - p)p- pp-^ + pP+ pP+\ (9) 

where the parameter p depends on the overlap of the two- 
photon amplitudes contributing to coincidence detection 
(r-r or t-t) and thus on the delay r. Ideally, p is zero if 
t = 0, increases as r increases, and approaches 1/2 when 
t — > 00. Figure 0] shows that this model largely explains 
the experimental results, except that 0\ is not actually 
equal to O2 as in the model. The reason is uncertain, 
but this inequality may be due to polarization-dependent 
1Z/T ignored in the model. 

We conclude that the function of the ideal BS is de- 
scribed by V~ , which depends on the p-s phase shift as 
well as the coefficients TZ and T. These parameters obvi- 
ously depend on the detailed design of the BS. Properly 
designed BS would satisfy 1Z/T = 1 and A — T = si- 
multaneously, at least for photons that have a particular 
wavelength and are incident at 45°. In this idealized case, 
V~ agrees with one of the four Bell states: a triplet state. 
This shows that the common belief that an ideal BS be- 
haves as a singlet-state filter is a misconception due to 
overlooking the helicity reversal of the photon upon re- 
flection at the BS. 

In summary, we used QPT to investigate the nonideal 
two-qubit quantum-state filter and identify its function. 
We showed that our results are satisfactorily explained 
by a simple model with several reasonable assumptions 
that a common belief about the function of a BS followed 
by post-selection process is mistaken. 

We thank S. Ishizaka and A. Tomita for their helpful 
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FIG. 3: The process matrices of the state filter obtained for delays r of (A) 0, (B) 100, and (C) 350 fs. Only the real parts of 
the matrices are shown (in arbitrary units). Vertical arrows indicate the diagonal elements (x F )i5i5- 




FIG. 4: Calculated process matrices corresponding to Fig. [3] (see text). The following parameters were used: (A) p — 0.14, 
(B) p = 0.325, (C) p = 0.5 and the common parameters TZ/T = 0.76, 9x « 0.41tt, and 8 2 « 0.076tt. 
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